Stiffness of the Edwards- Anderson Model in all Dimensions 
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A comprehensive description in all dimensions is provided for the scaling exponent y of low-energy 
excitations in the Ising spin glass introduced by Edwards and Anderson. A combination of extensive 
numerical as well as theoretical results suggest that its lower critical dimension is exactly di = 5/2. 
Such a result would be an essential feature of any complete model of low-temperature spin glass 
order and imposes a constraint that may help to distinguish between theories. 
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Imagining physical systems in non-integer dimensions, 
such as through an e-expansion near the upper Q or 
lower critical dimension |2|, has provided many impor- 
tant results for the understanding of the physics in real- 
istic dimensions. Often, peculiarities found in unphysical 
dimensions may impact real-world physics and enhance 
our understanding 3]. Here, we will explore the varia- 
tion in dimension of the low-temperature behavior in the 
Edwards- Anderson (EA) spin glass model [4J. 

A quantity of fundamental importance for the mod- 
eling of amorphous magnetic materials through spin 
glasses is the "stiffness" exponent y 0, • As Hooke's 
law describes the response in increasing elastic energy 
imparted to a system for increasing displacement L from 
its equilibrium position, the stiffness of a spin configu- 
ration describes the typical rise in magnetic energy AE 
due to an induced defect-interface of size L. But unlike 
uniform systems with a convex potential energy func- 
tion over its configuration space (say, a parabola for the 
sole variable in Hooke's law), an amorphous many-body 
system exhibits a function more reminiscent of a high- 
dimensional mountain landscape || . Any defect-induced 
displacement of size L in such a complicated energy land- 
scape may move a system through many ups-and-downs 
in energy AE. Averaging over many incarnations of such 
a system results in a typical energy scale 



<\AE\>~L V (L-m»). 



(1) 



The importance of this exponent for small excitations 
in disordered spin systems has been discussed in many 
contexts H H Q, ©, M El E3 

Spin systems with y > provide resistance ("stiff- 
ness") against the spontaneous formation of defects at 
sufficiently low temperatures T; an indication that a 
phase transition T c > to an ordered state exists. For in- 
stance, in an Ising ferromagnet, the energy AE is always 
proportional to the size of the interface, i. e. y = d — 1, 
consistent with the fact that T c > only when d > 1. 
When y < 0, a system is unstable (such as the d = 1 
ferromagnet) to spontaneous fluctuations which prolifer- 
ate, preventing any ordered state. Thus, determining the 
"lower critical dimension" di^wheie yd, = 0, is of signif- 
icant importance [6J, ll-'l , and understanding the 



mechanism leading to di, however un- natural, provides 
definite clues to the origin of order |2(. For instance, in 
homogeneous systems with a continuous symmetry |l6j . 
such as the Heisenberg ferromagnet, the possibility of 
"soft modes" (Goldstone bosons) perpendicular to the 
direction of magnetization have proved to weaken order 
further, since defect-interfaces can broaden over the en- 
tire length L of the system (spin wave) to reduce the 
energy by a factor ~ 1/L. This effect manifests itself in 
a reduced stiffness, y = d — 2, and an increased di = 2. 

The addition of quenched disorder drastically compli- 
cates the nature of the ordered state since fluctu- 
ating interfaces can attain highly irregular geometries to 
take advantage of heterogeneities. In this case, y, and 
possibly di, may become anomalous. For instance, it was 
found that adding even a small random field (of zero 
mean) to an ordinary Ising ferromagnet destabilizes or- 
der at least in d = 2, and considerable insight into the or- 
dered state was gained in proving that di = 2 [Lsl llflL Enj l . 

For spin glasses, even a numerical determination of y 
is highly nontrivial when d > 2 ;2jj, requiring repeated 
solution of some of the hardest-known combinatorial op- 
timization problems |2*3 |. Thus, simulations have been 
rather limited in system size L and yielded widely var ying 
results with non-overlapping errors, even for d = 3 23 1. 
The consensus from these studies is that 2 < di < 3 24 1 , 
where yd may behave in a discontinuous manner 25 1. 
such that di could be integer as is the case for other sta- 
tistical models. On the other extreme, above the "upper 
critical dimension," well-known to be d u — 6 0, where 
fluctuations become irrelevant, numerical results can be 
compared with analytical predictions for critical expo- 
nents from mean-field theory |2fij . 

In this Letter, we consider an analytic continuation of 
V = Vd to real dimensions d > 0, based on extensive nu- 
merical simulations for y d in integer dimensions [iji |2lJ • 
That data suggest that yd is varying smoothly across di- 
mensions and that a straightforward low-order fit would 
provide a very accurate description. The validity and ac- 
curacy of this fit is gauged by how well it extrapolates 
to an exactly known feature, e. g. y\ = — 1. The fit has 
a zero at di rs 2.4986 and yields y 5 / 2 ~ 0.0008; strong 
evidence that di = 5/2. This value can be corroborated 



2 



by a similar but independent fit of the existing data for 



T g |28|, |29j, |30j . A speculative calculation based on mean- 
field arguments [3l[ , recently put on a more rigorous ba- 
sis [3^, has predicted the disappearance of replica sym- 
metry breaking (RSB) j3j|, taken to be a characteristic of 
low-T glassy order, for dimensions d < di — 5/2 exactly. 
Our results here lend significant support to the low-T 
structure of spin glasses advanced in Refs. Isil ls^ . 

Applying a new optimization heuristic |34j to spin 
glasses on bond-diluted lattices [23j with up to 10 6 spins 
allowed us to present values for yd in d — 3, . . . , 7 that 
are improved or are entirely unprecedented: 2/3 = 0.24(1), 
2/4 = 0.61(1), 2/5 = 0.88(5), y 6 = 1.1(1), and 2/7 = 1-24(5). 
The precision of these values originates from extensive 
numerical simulations [2]} over a large number of sys- 
tems sizes and bond fractions, yielding 30 or more data 
points significant for scaling in each dimension, with each 
point resulting from an average over 10 4 — 10 6 instances. 

We can supplement this list of values for yd by those 
for d = 1 and 2, assuming a bond distribution that 
in addition to being symmetric and of unit width also 
has to be continuous at the origin. Only for yd > 
is its value independent of details of the bond distribu- 
tion , where the energy scale in Eq. Q diverges. For 
d < 2, discrete ± J bonds appear to lead to trivial scal- 
ing, yd<2 = 0, H3- whereas continuous bonds gen- 
erally provide nontrivial results [24j |. For instance, in a 
d = 1 spin chain of length L, any defect at T — af- 
fects only the absolute weakest bond. Hence, for ±J 
bonds, or any other distribution bounded away from the 
origin, the absolute weakest bond is a non-zero constant 
for large L. For a distribution continuous at the origin, 
the absolute smallest of L bonds scales as ~ thus 
2/1 = —1 exactly. The value for 2/9 has been considered 
repeatedly over the years 0, @, yJj |36|, |33 , and has been 
determined with some consistency to be 2/2 = —0.282(2) 
for continuous bonds. The accuracy for 2/2 is expected to 
be significantly better than for any yd with d > 2, since 
finding exact ground states for d = 2 spin glasses is ac- 
complished with polynomial time algorithms, facilitating 
larger sizes and better statistics. 

From these values, a clear picture for the stiffness ex- 
ponent emerges. As Fig. Q] shows, the exponent appears 
to vary smoothly with dimension d, at least between 

1 < d < 6. (It can be expected that yd would develop a 
cusp at the upper critical dimension d u — 6.) Hence, we 
expect that a simple polynomial fit of the numerical data 
may provide an accurate interpolation for yd- In light of 
the much reduced uncertainty in the value of 2/3 com- 
pared to earlier results, which ranged between 0.19(1) to 
0.27 H3, such a fit should provide an accurate prediction 
for di, i. e. the zero of yd- To fit the data, we proceed 
in the following way: We fit only the numerical data for 

2 < d < 6, weighted by their error bars; the quality of the 
fit is judged on the basis of how well the exactly known 
data point y\ = — 1 is reproduced. We find that the opti- 
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FIG. 1: Plot of the spin glass stiffness exponent yd as a func- 
tion of dimension d. Shown are the data for hypercubic lat- 
tices from Refs. [23ll35jl . a weighted cubic fit to the yd data on 
2 < d < 6, and the RSB prediction from Ref. [H in Eq. ©. 
The fit predicts to 0.1% that the lower critical dimension de- 
fined through yd l — is di —5/2 (vertical line). It reproduces 
the exact result, y\ = — 1 (checkered square), to 0.8% and sug- 
gests yo — —2. Above the upper critical dimension, d u = 6, 
the data significantly deviates from Eq. (JSJ (solid line). 



mal fit is of cubic order: a square fit does not even have 
enough variability to match all the included data well, 
let alone 2/1; fits of higher order become unstable since 
the number of fitting parameters approaches the number 
of data points. The resulting cubic fit, 



Vd 



-1.988 + 1.125(2 - 0.1533d 2 + 0.0086d 3 , (2) 



has a number of desirable properties. In particular, the 
coefficients for increasing powers of d rapidly decrease, 
testifying to the stability of the fit. Consequently, Eq. @ 
predicts the exactly-known result, 2/1 = —1, to less than 
0.8% and further extrapolates to 2/0 = —2 (within 0.7%). 
The fit has a zero at d t « 2.4986 and yields y 5/2 » 0.0008; 
suggesting that di = 5/2. 

In fact, a replica theory calculation on a lattice in ar- 
bitrary dimensions [3l| yields di =5/2 exactly. It uses 
a Landau-like expansion in the neighborhood of T g for 
the order parameter function for two coupled replicas. 
Both replicas are held at an identical equilibrium state 
on one open boundary, and held at two distinct states on 
the opposite boundary. While the "defect" thus created 
leads to a different exponent, say, ad ^ yd, their zero's 
should coincide. The defect free energy AF signals that 
the distinction between overlaps becomes irrelevant ex- 
actly in d — 5/2, where replica symmetry gets reestab- 
lished throughout. The ad- hoc free-energy extremization 
in Ref. 13 ll has recently been put on a more rigorous 
basis 32], using the inverse range as an expansion pa- 
rameter near a saddle point of a lattice spin glass with 
variable interaction range. It is argued that this AF be- 
haves similar to that in ordinary ordered systems with a 
continuous symmetry, such as the Heisenberg ferromag- 



3 



net, in which the Goldstone mode undermines order al- 
ready in di = 2 0| . Yet, in the Ising spin glass order is 
further weakened by the bond disorder without any fer- 
romagnetic bias, making any contribution to the defect 
energy sub-additive ((AF) = but (|AF|) > Oh hence 
ad = d— 2— |. The replica calculation in Refs. |3J,|33| is 
presently the only realistic theory with an analytic pre- 
diction for di. Considering the simplicity of a di = 5/2, 
one may wonder how essential the replica argument is. 

We can put these predictions further into context 
by compa ring with the Migdal-Kadanoff approximation 
(MK) |33, |39j , which has proved to resemble many fea- 
tures of low-dimensional spin glass quite well |40j. In 
MK, a d-dimensional lattice is represented by a hierar- 
chical graph where each bond of range sx L between two 
spins consists recursively of 6 parallel series of s bonds of 
range L, leading to 0, Oil 
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In 6 
In s 



(3) 



Recent simulations of MK in non-integer d (using 6 = 3 
and continuous s) have come surprisingly close to the 
same result 0, di ~ 2.5. Already Ref. [f| found analyt- 
ically for MK, within a Gaussian approximation, that 



,.MK 
Vd 



d-l In (1 



21n2 



(4) 



which has been improved to a d — > 00 expansion of yd in 
Ref. 0. Clearly, Eq. (0J fails at d = 1, and its linear 
form misses the essential features of the data in Fig. ^ 
but it provides a decent estimate for d = 2 and 3, and pre- 
dicts df IK w 2.46. (Similarly, a suggestion from Ref. [Ilj] 
that [d — l)/2 — y = const is ruled out by the data.) 

Yet, MK predictions are generally ambiguous, since 
various combinations of 6 and s in Eq. J3J can represent 
the same d. For instance, for s> 1 each series of bonds 
at T = can be replaced by the (absolute) smallest bond 
(like the d = 1 spin chain above). Drawn from a Gaussian 
distribution of zero mean and unit variance, the equiva- 
lent bond replacing each series is of typical size ~ 1/s and 
has a random sign. Putting b such iid bonds in parallel 
thus yields another effective bond of Gaussian distribu- 
tion with zero mean and width ~ y/b/s. Hence, for an 
interfacial energy AE to be scale invariant [i. e. y = in 
Eq. 0J] requires b = s 2 for b, s — ► 00, implying df IK = 3 
according to Eq. 0. This result is at odds with compu- 
tations for 6 = 3 |24j and variable s. Furthermore, using 
s = 2 and an analytic continuation for 6 — > 1/2, one can 
show that y^ K = — 1; quite inconsistent with Fig.^ 

As an independent test for the picture emerging from 
Fig. ^ we study the values for the glass-transition tem- 
peratures as a function of dimension, T g (d), for d > 3. 
To this end, we have gathered the latest values for T g (d) 
(usin g ± J bonds) from the literature, based on simula- 
tions 1H, E3, or high-T series expansions |30l . We 




FIG. 2: Plot of the glass transition temperature T g (d) as a 
function of dimension d. Data points are drawn from various 
studies IH, H||3(3, ^ The solid line represents the fit 
in Eq. ©. That fit captures the known asymptotic behavior 
T g (d) ~ V2d for d — > 00 and predicts a zero at d ss 2.491. 



found T g {3) w 1.15, T 3 (4) » 2.03, T 9 (5) » 2.548, T g {6) » 
3.026, T g (7) w 3.385, T s (8) w 3.694, and T g (9) w 4.011. 
Determining di by fitting T g (d) is somewhat more com- 
plicated, since it vanishes for all d < di (unlike yd) and 
Ref. [3 argues that T g {d) - y/d-di for d ^ dj . On the 
other hand, for the Sherrington-Kirkpatrick (SK) model 
T^ K = J SK 0| and J ~ J SK /V2d, so T ff (d) - V2d for 
d — > 00 and J = 1. Satisfying these constraints, we fit 



T ff (d) pa v/2(d- 2.491) [ 1 + 



1.237 2.425 



d 2 



(5) 



which is plotted in Fig. [21 Again, incorporating more 
terms in the fit leads to drastic oscillations in the coeffi- 
cients. While the predicted result, di sa 2.491, is not quit e 
as accurate or stable as above, it is quite consistent. |51| 
Unlike the consistency between mean-field arguments 
and numerics for di, we can observe a large discrepancy 
for d > d u = 6. An RSB calculations from Ref. 
around the mean-field limit (rf — > 00) 45] predicts 



y d = d(l-p) fa - 



(d > d u = 6), 



(6) 



where the exponent p, predicted [ill l26j | and mea- 
sured [3, EE E3 to be « 3/4, describes the width of the 
distribution of ground state energies in the SK model. As 
Fig- n shows, the finite-dimensional results substantially 
disagree on this point with replica theory. A plausible 
interpolation may be yd ~ (d — 2)/4, which approaches 
Eq. JSJ for d — > 00. 

In conclusion, we have provided a series of independent, 
mutually consistent arguments for a lower critical dimen- 
sion of di = 5/2 for the Edwards- Anderson spin glass. 
This self-consistency has a number of important implica- 
tions. For one, it validates the numerical values for yd for 
integer dimensions as quoted above [2E H3 ; any change 
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beyond the error bars of 1/3 in particular would signif- 
icantly distort the obtained di. The accuracy of those 
values, in particular y e and j/7, makes the comparison 
with the mean-field prediction for d > 6 relevant, al- 
lowing for a rare direct comparison between lattices and 
mean- field predictions across the dimensions. In turn, 
the apparent square-root singularity near di validates the 
arguments from Ref. for the behavior of T g . Most 
importantly, the consistency of these numerical results 
with the analytical treatment in Refs. lend cred- 

ibility to its assumptions and implications. That theory 
assumes replica symmetry breaking to exist below T g for 
finite-dimensional lattices, which seems incompatible 48] 
with the droplet scaling picture developed in Refs. 0,Ta|- 
It further derives that the treatment is unaffected by 
the existence of an external field, hence predicting an dc 
Almeida-Thouless line, in contrast with some recent ex- 
perimental 0| and numerical findings [H(J • Our findings 
here suggest a careful reexamination of those ideas. 
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